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^ Abstract 

The present note deals with a nonstandard systems of differential equations de- 
^ . scribing a two-species phase segregation. This system naturally arises in the asymp- 

. totic analysis carried out recently by the same authors, as the diffusion coefficient 

in the equation governing the evolution of the order parameter tends to zero. In 
particular, an existence result has been proved for the limit system in a very general 
framework. On the contrary, uniqueness was shown by assuming a constant mobil- 
ity coefficient. Here, we generalize this result and prove a continuous dependence 
property in the case that the mobility coefficient suitably depends on the chemical 
potential. 
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Continuous dependence for a nonstandard phase field system 



1 Introduction 

In this paper, we address the system 



(1 + 2g{p)) dtfi + fig'{p) dtp - div(K(/i)V/i) = (1.1) 

dtp + f'ip) = fig'ip) (1.2) 

{k{p)V p) ■ u\r = (1.3) 

/x(0) = po and p(0) = po (1.4) 



of differential equations and boundary and initial conditions in terms of the unknown 
fields p and p; equations f ll.ip -f lTT^ are meant to hold in a bounded domain C with 
a smooth boundary F and in some time interval (0, T), and u in fll.3p denotes the outward 
unit normal vector to F. The recent paper [2] investigated the existence of solutions to 
the above system: actually, a solution was found by considering the analogous system 
in which the ordinary differential equation (11. 2p is replaced by the partial differential 
equation 

dtp — crAp + /'(p) = pg\p) with the boundary condition dup\T = 0, (1.5) 

and then performing the asymptotic analysis as the diffusive coefficient a tends to zero. 

Such a modified system arises from the model introduced in [5], which describes the 
phase segregation of two species (atoms and vacancies, say) on a lattice in the presence 
of diffusion. It turns out to be a modification of the well-known Cahn-Hilliard equations 
(see, e.g., [3 [8]). The state variables are the order parameter p (volume density of one of 
the two species), which must of course attain values in the domain of the nonlinearities 
g' and /', and the chemical potential p, which is required to be nonnegative for physical 
reasons. The initial-boundary value problem for the PDE system has been studied in a 
series of papers with a number of obtained results: here, we confine ourselves to quote 
the former [2 El [6] and latter [21 [3l [1]. 

In the mentioned papers, the function g is taken as a smooth nonnegative and possibly 
concave function (like it looks here), while the function / represents a multi-well potential: 
in this respect, a thermodynamically relevant example for / is the so-called logarithmic 
potential, in which /' is given by the formula 

/'(p)=lni±^-2cp for pG (-1,1), (1.6) 
1 - p 

with c > 1 in order that / actually presents a double well. The class of the admissi- 
ble potentials may be rather wide and include both the standard double-well potential 
defined by 

/(p) = i (p2 - 1)2 for PGM (1.7) 

and potentials whose convex part /i is just a proper and lower semicontinuous function, 
thus possibly non-differentiable in its effective domain. In such a case, the monotone 
part f[ of /' is replaced by the (possibly) multivalued sub differential dfi and (11. 5p has 
to be read as a differential inclusion. In [2] , this wide class of potentials was considered. 
Moreover, in [2] the mobility coefficient n in (11. ip and (11. 3p was allowed to depend also 
on p. 



Colli 



GiLARDI 



Krejci 



Sprekels 



3 



Therefore, the existence result for system fll.ip - fll.4p proved in [2] is very general. 
On the other hand, the solution constructed in this way is rather irregular, in principle 
(due to a lack of regularity for /i). Nevertheless, it has been shown to be unique (and a little 
smoother than expected) provided that the mobility coefficient k is a positive constant. 

The aim of the present paper is to generalize the uniqueness proof performed in |2] to 
the case of a mobility coefficient depending on the chemical potential, exactly as in (11. ip 
and (II. 3p . Moreover, the continuous dependence of the solution on the initial data is 
shown in terms of suitable norms. Of course, in order to accomplish our program, a 
natural uniform parabolicity condition is required for k. 

The paper is organized as follows. In the next section, we list our assumptions and 
rewrite problem (Il.ip -( |L4|) in a precise form. In Section |3l we state and prove our 
uniqueness and continuous dependence result. 



2 Assumptions and notations 

We first introduce precise assumptions on the data for the mathematical problem under 
investigation. We assume f2 to be a bounded connected open set in with smooth 
boundary F and let T G (0, +oo) stand for a final time. We set for brevity 



The symbol ( ■ , ■ ) denotes the duality product between V*, the dual space of V, and V 
itself. For the nonlinearities we assume that there exist real constants k*, p*, p*, C,*, 
and ^* such that the combined conditions listed below hold. 



Notice that important potentials like (II. 6p and (I1.7P fit the above requirements with 
suitable choices of g and of the constants. For the initial data, we require that 



V:=H\n), H:=L\n), and g:=fix(0,T). 



(2.1) 




(2.2) 
(2.3) 
(2.4) 
(2.5) 
(2.6) 
(2.7) 
(2.8) 
(2.9) 
(2.10) 



fio eV n L°°{Q) and po > a.e. in Q 
Po ^ V and p* < Po < p* a.e. in Q. 



(2.11) 
(2.12) 



Now, we recall the part that follows from the asymptotic analysis performed in [2] and 
is of interest for the present paper. 
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Theorem 2.1. Assume that (12.21) - fl2.12p hold. Then there exists at least one triplet 
(/i, p, C.) that satisfies 

fxe L°°iO,T;H)r]L^{0,T;V)nL°°{Q) and fi > a.e.mQ (2.13) 

p G L°°(0,T;\/), dtpeL°°{Q), and < p < p* a.e.mQ (2.14) 

eGL°^(g), ee/3(p) and < ^ < T a.e.mQ (2.15) 

M := (1 + 2^(p))p G Vr^'P(0,T;r*) nL2(0,r;iy^''?(fi)) (2.16) 

/or some p, g > 1 and solves the problem 

{dtu{t),v)+ I K{p{t))Vfi{t)-Vv= [ p{t)g'{p{t))dtp{t)v 
Jn Jn 

for allv eV and a.a. t e {0,T) (2.17) 
9tP + ^ + T^ip) = fJ- g'ip) a.e.inQ (2.18) 
u(0) = (l + 2(7(po))po and p(0) = po a.e.inVt. (2-19) 

Remark 2.2. We notice that fl2.17p actually is a weak form of equation (11. ip (with the 
boundary condition (11.31) since the test function v is free on the boundary). Indeed, 
whenever p is smoother with respect to time, one can compute dtu by the Leibniz rule 
and see that the differential equation hidden in the variational equation (I2.17P coincides 
with (II. ip . We also observe that [2] precisely yields p = 4/3 and q = 3/2 in (12.161) . On 
the other hand, the regularity u G L'^{0, T; V) follows immediately from u = (l + 2g[p)^ji 
thanks to (I2.13p - (l2.14p . whence one can take q = 2. 

In [2], an additional result was proved that deals with continuous dependence on the 
initial datum po. Here, we adapt the statement to our purposes. Indeed, we also consider 
possibly different initial data for the chemical potential (even though they do not enter 
the final estimate, directly). 

Proposition 2.3. Assume that (I2.2l) - (l2.12p hold. Let {po,i,Po.i), ^ = 1,2, be two sets of 
initial data satisfying (I2.1ip - (I2.12I) . and let {pi,pi,C,i), i = 1,2, be two solutions to the 
corresponding problem ( 12.17P -( !2.19p that satisfy the regularity assumptions ( ]2.13p -( !2.16p . 
Then the following estimate holds true: 

\{P1 - P2)(t)| + l{\9tiPl - P2)\ + 1^1 - 6l) (S) ds 

< C (^\po,i - PoM + ^ (l/ii - P2I + (1 + Pi)\pi - P2I) (s) ds^ (2.20) 

for every t G [0, T] and a.e. in Q, where C depends only on the constants and the functions 
mentioned in our assumptions (12. 2p - (12.100 on the structure of the system. 



Proof. In order to give just an idea how to obtain (12.201) . let us point out that the pro- 
cedure consists in testing the difference of two equations (12.181) by sign(,^i — ,^2)- Indeed, 
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setting Wi = Hig'{pi) — TT{pi), i = 1,2, and multiplying the identity 

dtipi - P2) + i^i - ^2) = wi - W2 (2.21) 
by sign(,^i — .^2), it is not difficult to infer that 

dt\pi - P2\ + - ^2\ < {wi - W2\ a.e. in Q. (2.22) 



r 16 -61(5) ds 
Jo 



Thanks to the Lipschitz continuity properties in (12. 8p . and integrating (I2.22p only with 
respect to time, we obtain that for t G (0, T) it holds 

\pi - P2\it) 

Jo 

< c ^|po,i - Po,2\ + (l/ii - fJ'2\ + (1 + Ati)|Pi - P2\){s)ds 
a.e. in Q. Moreover, note that (12.211) implies 

/ dt\pi - p2\{s) ds < / {\wi - W2\ + - ^2\){s) ds, 
Jo Jo 

whence (I2.2UI) easily follows. □ 

In [2] , the uniqueness of the solution given by Theorem 12.11 (as well as the regularity 
dtfi G L'^{Q)) was proved under an additional assumption, namely: 

Theorem 2.4. Assume (12.21) - (12. 12p and that k is a positive constant. Then the solution 
{fj,,p,^) given by Theorem \2.1\ is unique. 

The aim of this paper is to improve this result by showing that uniqueness and con- 
tinuous dependence hold in the more general framework of Theorem 12. H as stated in the 
forthcoming Theorem 13. 1[ 



3 Uniqueness and continuous dependence 

In this section, we prove the uniqueness and continuous dependence result for the solution 
to problem fl2Tri) -f l239|) stated below. 

Theorem 3.1. Assume that the conditions (I2.2p - (l2.12p are satisfied. Then the solution 
{p,p,^) given by Theorem \2.1\ is unique. Moreover, let {po,i,Po,i), i = 1,2, be two sets 
of initial data satisfying (I2.11l) -( l2.12p . and let {pi,pi,C,i), i = 1,2, be the corresponding 
solutions, which fulfill (12.170 -( l2.19p with Pq = po,i o,nd po = po,i> i = 1,2. Then there 
exists a constant C , depending on the data through the structural assumptions, such that 

IIPi — P2||L2(0,T;/f) + IIPi " P2\\l°°(0,T;H) + ||6 " ^lU^Q) 

< C {||/io,i - P'oAh + ||po,i - PoAh} ■ (3.1) 
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Proof. We just prove continuous dependence, with uniqueness as a byproduct. Through- 
out the proof, we account for the well-known Holder inequality and for the elementary 
Young inequality 

ab < ea^ + —b"^ for every a,b > and e > 0. 

Moreover, in order to simplify the notation, we use the same symbol small-case c for 
different constants, which may only depend on Q, the final time T, the nonlinearities k, 
/, g, and the solutions under consideration. Thus, the meaning of c may change from 
line to line and even within the same chain of inequalities. In contrast, we choose capital 
letters to denote precise constants we want to refer to. Finally, we set 

Qt:=nx (0, t) for t E (0, T]. (3.2) 

Our argument relies on a suitable adaptation of the technique developed in [2] with the 
help of the function : [0, +oo) — )■ M defined by 

K{m) := / K{m') dm' for m > 0. (3.3) 
Jo 

We have indeed K' = k, whence VK{fi) = K(/i)V/i, so that (I2.17P becomes, with k : = 



{dtu{t),v)+ / Vk{t)-Vv= / fi{t)g'{p{t))dtp{t)v 
Jn Jn 

for allv eV and a.a. t e (0, T). (3.4) 

We remark at once that (12.31) yields that for every mi , m2 > it holds 

n^\mi — < \K{mi) — K{'m2)\ < K,*\mi — (3-5) 
{nil - "^2) {K{mi) - K{m2)) > K,*{mi - m2)^. (3.6) 

In our proof, we use the equation obtained by integrating ( 13. 4p with respect to time rather 
than (13. 4p itself; namely, for every v eV and t E [0,T] we have 

u{t)v- [ u{0)v+ [ Vk{t)-Vv= [ (f;^fi{s)g'{p{s))dtp{s)ds)v, (3.7) 
n Jn Jn Jn^ ^ 

where ^ ^ 

k{t):= [ k{s)ds= [ K{fx{s))ds, (3.8) 
Jo Jo 

and where n(0) = (1 + 2g{pQ))pQ according to the first Cauchy condition (I2.19p . It is 
worth observing that the pointwise values of u in the integrals over Q are well defined. 
Indeed, the boundedness of u (derived from the boundedness of fi and g{p)) and the 
regularity of dtu (cf. f l2.13p -( l2.14p and f l2.16p ) ensure that u is weakly continuous, e.g., as 
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an if-valued function. Now, we let (/io.ii Po,j); = l!2, be the initial data and pick two 
solutions {ni, Pi, ^i), i = 1, 2. Then, let us define the corresponding functions Ui, ki, and ki 
(according to fl2.16p . (13 ■4p . and (I3.8P ). as well as the new ones ■ji, as follows: 

7i := 1 + 2g{pi) , Ui := 7i fii , ki := K{pi) , 

k^{t) := / A;i(s)cis fortG[0,T], z = l,2, 

in order to simplify the notation. For the same reason, we set 

Po ■= Po,i ~ Po,2 ; Po •= Po,i ~ Po,2 
p:=pi-p2, p:=Pi-p2, ^:=6-6 

7 := 7i — 72 , u := Ui — U2 , k := ki — k2 , and k := ki — k2 . 

At this point, we write (13. 7p at the time t for both solutions and choose v = k{t) = dtk{t) 
in the difference. We obtain 

/ u{t)k{t)+ I Vk{t) -Vdtkit) 
Jq Jn 

= I ((l + 2(?(po,i))po + 2po,2(^7(po,i)-^7(Po,2)))A;(t) 
Jn 

+ j^[f^{pl{s)g\p^{s))^Ms) - P2{s)g\p2{s))dtP2{s))dsyit). (3.9) 

We estimate each term of separately. By accounting for (Q-d^SD, (127[3D - (I27[4D . 
as well as for (I3.5p - (I3.6I) . we have 

uk = (71P1 — 72^2)^ = 7iP^ + 7P2^ > /t*|p|^ — c\p\ \p\ a.e. in Q, whence 

/ u{t)k{t)>K^ I \p{t)\'-C I |p(t)||/i(t)|>^ / |p(t)r-C / \p{t)\\ 

Jn Jn Jn ^ Jn Jn 

Next, we clearly see that 

f vk{t).vdMt) = lj^ I m{t)\\ 

With the help of ([2l]), fl2TT]) - (l212D and ([33D, we can control the first term on the 
right-hand side of (13. 9p : 

/ ((1 + 2(7(po,i))po + 2/io,2((7(po,i) - g{p(^,2)))k{t) 
Jn 

<c [ (|po| + |po|)|p(t)|<c(||po|lH+IIPo||H) + T / l/^^^)!'- 
Jn ^ Jn 

In order to estimate the second term, we observe that 
\P'ig'ipi)dtpi - P2g'ip2)dtp2\ 

< \P\ \g'{pi)\ \dtpi\ + IP2I l^'(pi) - ^'(P2)| \dtpi\ + IP2I |^'(P2)| \dtp\ 

< c{\p\ + \p\ + \dtp\) a.e. in Q, 



8 



Continuous dependence for a nonstandard phase field system 



thanks to our regularity assumptions on the solutions and on the structure (cf. fl2.13p - 
(I2.14P and (12. Sp ). By owing to Proposition I2.3[ we deduce that 
ft 

{l^iis)g'ipiis))dtpi{s) - H2{s)g\p2{s))dtp2{,s)) ds 



< c {\p{s)\ + \p{s)\ + \dtp{s)\) ds < c[\po\ + £{Hs)\ + \p{s)\) dsj a.e. in fi. 
Therefore, we have 

Io{f^i{s)9'{pi{s))dtPi{s) - fi2{s)g'{p2{s))dtP2is)) ds)k{t) 



ij<" 

<c f (\Po\+ IoM')\ + \p{s)\) ds) \p{t)\ 



< 



^ / \m' + c[ {\Po\'+ IoHs)\'ds + J^\p{s)\'ds} 

Ht)\' + c\\pofH + C [ \p\' + c[ \p\'. 



T 



Qt 



By combining the above equalities and inequalities with fl3.9p . we infer that 

^* f I / . \ 1 9 1 

T 

and an integration with respect to time yields 

f^7k / I lO 1 



iM^)r+^^/ ivMt)r<c(ii/ioii?,+iiPoiii)+c/ \p{t)\'+c[ i/ir+c/ ipp 



JQt ^ 



4 

<c(||/io||?,+ ||po||?,)+c/ \p\'+f([ l/iP+Z IpP)^^ 

JQt ^0 ^JQs JQs ' 

<c(||/io||^+||po||i)+c /V/ \p\')ds + cf \p\'. (3.10) 

^0 ^^Q, ^ JQt 

Now, let us consider (I2.20p . Squaring and integrating over fl, then applying Holder's 
inequality on the right-hand side, we easily obtain that 



\p{t)\' < D\\po\\l + D / l/ip + D / \p\' (3.11) 

JQt JQt 

for some positive constant D. Moreover, by integrating (I2.20p over Q and then squaring, 
we arrive at ^ 

( [ \^\) <D\\p4j, + D I \p\' + D f \p\\ (3.12) 
\JQt / JQt JQt 

where D is the same constant as before, without loss of generality. Hence, we multiply 

(I3.10p by 12£)/k* and add it to (13. lip and (I3.12p . This computation leads to 

DMiHQt) + Mt)rH + m\UQ.) 

<c(||/io|li+||po||?,)+c / Ml^Q^)ds + c [ \\p{s)\\lds. (3.13) 
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At this point, it suffices to apply the Gronwall lemma to deduce a variation of ( 13.11) with 
the squared norms. Therefore, (13. ip is completely proved. □ 

Remark 3.2. Clearly, just a few of the assumptions f l2.2p -( l2TT2l) are used in the above 
proof. The whole set of hypotheses has been listed in the statement of Theorem 13.11 in 
order to ensure both the existence of a solution satisfying fl2.13p - fl2.16p and the validity 
of estimate fl2.20p . according to Theorem 12. II and Proposition 12.31 
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